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Abstract 

Experimental data on pressure distribution and 
heat transfer on a turbine airfoil were obtained over a 
range of Reynolds numbers from 0.75 to 7.0x10 and a 
range of turbulence intensities from 1.8 to about 15 per- 
cent. The purpose of this study was to obtain fundamen- 
tal heat transfer and pressure distribution data over a 
wide range of high Reynolds numbers and to extend the 
heat transfer data base to include the range of Reynolds 
numbers encountered in the Space Shuttle main engine 
(SSME) turbopump turbines. Specifically, the study 
aimed to determine (1) the effect of Reynolds number on 
heat transfer, (2) the effect of upstream turbulence on 
heat transfer and pressure distribution, and (3) the rela- 
tionship between heat transfer at high Reynolds num- 
bers and the current data base. The results of this study 
indicated that Reynolds number and turbulence intens- 
ity have a large effect on both the transition from lami- 
nar to turbulent flow and the resulting heat transfer. 
For a given turbulence intensity, heat transfer for all 
Reynolds numbers at the leading edge can be correlated 
with the Frossling number developed for lower Reynolds 
numbers. For a given turbulence intensity, heat transfer 
for the airfoil surfaces downstream of the leading edge 
can be approximately correlated with a dimensionless 
parameter. Comparison of the experimental results 
were also made with a numerical solution from a two- 
dimensional Navier-Stokes code. 

Nomenclature 

A area, ft 2 

B curve-fitting polynomial coefficient 

C chord length, ft 

d leading edge diameter, ft 

E mean voltage, volts 

Fr(0) Frossling number at stagnation point, 

Nu d /(Re inlet ) a5 

h heat transfer coefficient, Btu/hr/ft 2 /°F 


k thermal conductivity of air, 

«0.0147 Btu/(hr/ft/°F) 

L total length of airfoil pressure or suction 

surface, in. 

Nu Nusselt number, h d/k 

P pressure, psia 

Pr Prandtl number 

Q heat flow rate, Btu/hr 

Re ex it Reynolds number, based on axial chord and 

airfoil exit conditions 

R e iniet Reynolds number, based on leading edge 
diameter and inlet conditions 

r recovery factor, (Pr) 1 / 3 « 0.89 

5 blade pitch, in. 

s distance along airfoil pressure or suction 

surface from airfoil stagnation point, in. 

T temperature, °R 

Tu turbulence intensity, percent 

V fluid velocity, ft/sec 

W flow rate, lb/sec 

6 flow angle, rad or deg 

A x longitudinal integral length scale 

fi dynamic viscosity, lb/ft -sec 

p fluid density, lb/ft 3 

Subscripts: 

0,1, 2, 3, 4 curve-fitting polynomial coefficients 

c chord length 

d leading edge diameter 

g gas 

i input 


* 


Resident Research Associate, NASA Lewis Research Center. 


1 


in inlet 

L loss 

LC liquid crystal 

R recovery 

s surface distance from airfoil stagnation point 

st static 

t total 

Superscript: 

fluctuating velocity component 
Introduction 

The turbines used in the Space Shuttle main engine 
(SSME) run at very high pressures and Reynolds num- 
bers. At these extreme conditions, the heat transfer rates 
to the turbine airfoils are expected to be significantly 
higher than those of current aeropropulsion turbines. A 
detailed knowledge of heat transfer at these extreme con- 
ditions is necessary to predict airfoil surface tempera- 
tures and satisfy life goals. 

A large body of flow and heat transfer data over 
airfoils for current aeropropulsion turbines at lower 
Reynolds numbers on the order of 10 6 is available in the 
literature. However, very little information is available 
for the range of higher Reynolds numbers found in the 
SSME turbopump turbines. 

A Variable Reynolds Number Heat Transfer 
Cascade Facility was constructed at the NASA Lewis 
Research Center. Its purpose was to conduct a funda- 
mental study on heat transfer at Reynolds numbers 
ranging from current aeropropulsion gas turbine levels 
to those found in the SSME turbopump turbines. The 
effect of free-stream turbulence on heat transfer and the 
relationship between heat transfer data at high Reynolds 
numbers and the current data base was investigated. 
Turbulence-generating grids were used to vary the tur- 
bulence intensity levels in the test section. The facility 
is capable of operating over a range of Reynolds num- 
bers from 0.75 to over 7.0xl0 6 (based on axial chord 
and airfoil exit conditions) and a range of turbulence 
intensities from 1.8 to about 15 percent. 

Heat transfer and airfoil pressure distribution data 
were obtained at Reynolds numbers of 0.75, 1.5, 3.0, 5.0, 
and 7.0xl0 6 . The measured turbulence intensity levels 
ranged from 1.8 percent at a Reynolds number of 
7.0xl0 6 with no grid (clear tunnel case) to 15.1 percent 
at a Reynolds number of 0.75 XlO 6 with 1/4- in. grid. 
Preliminary results from this facility were reported in 
Schobeiri, McFarland, and Yeh. 1 This paper presents the 
latest heat transfer and pressure distribution results 


from this facility. Experimental heat transfer data are 
also compared with a numerical solution from the 
TRAFC2D 2 code. 

F acility 

General Description 

Figure 1 is an illustration of the test facility, 
including a cutaway view of the test section. Com- 
pressed air at 40 psig and ambient temperature flows 
through an 18-in.-diameter pipe which houses screens 
and perforated plates for flow conditioning. The maxi- 
mum flow rate of the facility is about 38 lb/sec; the air 
flow exits to the laboratory exhaust system at 1.93 psia. 
Upstream of the test section, boundary layer air is bled 
from four sides to provide a uniform velocity profile at 
the test section entrance. The main air, flow as well as 
each individual bleed air flow, is measured with a sharp- 
edged orifice. 

Test Section 

The test section with the cover plate removed is 
shown in Fig. 2. Visible in this photograph are the 
airfoil, the contoured sidewalls, the hot film anemometer 
probe, the inlet and exit flow angle measurement probes, 
and the Pitot static probe. A turbulence grid, not 
visible, is located behind the grid cover. Two turbulence 
grids were used to vary the turbulence intensity in the 
test section. The grids were of the square-mesh, bi-plane 
type with square bars, one having l/8-in. square bars 
and the other having 1/4- in. square bars. Tests were 
also conducted using a clear tunnel only (no grid). The 
ratio of mesh-to-bar width spacing was 4.5, giving a 
60-percent open flow area. The traversing anemometer 
probe, mounted on the test section, is located 3.6 lead- 
ing edge diameters (6.19 in.) downstream from the grid 
and 0.77 diameter (1.31 in.) in front of the airfoil 
leading edge. The anemometer probe was used to survey 
the turbulence intensity in the pitchwise direction. 

Airfoil 

The airfoil used in these tests has an axial chord of 
8.5 in., an actual chord of 11.09 in., a span of 8.5 in., a 
solidity of 1.32, a camber angle of 51.7°, and a two- 
dimensional, constant-section vane to provide a uniform 
flow field. Coordinates of the airfoil and other geometric 
information are given in Schobeiri, McFarland, and 
Yeh. 1 Two airfoils were fabricated: one was used for 
pressure distribution tests and the other for heat transfer 
tests. The pressure distribution airfoil is shown in the 
photograph of Fig. 3. Of the 90 pressure taps on the air- 
foil, 52 were located at the midspan, and 19 each were 
located at the hub and tip sections to record any three- 
dimensionality of the flow field around the airfoil. 
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The airfoil used for heat transfer tests is shown in 
Fig. 4. It is overlayed with a thin (0.001-in.) Inconel 
sheet. Passing an electric current through the Inconel 
sheet gave a constant heat flux boundary condition. The 
Inconel sheet was first coated with black paint, then 
sprayed with a liquid crystal and clear binder mixture. 
Grid lines were drawn on the airfoil surface to locate the 
temperature field in the data reduction process. White 
dots were added as an aid to locate the temperatures 
indicated by the liquid crystals. Details on the use of the 
liquid crystal technique on airfoils are given by 
Hippensteele, Russell, and Torres. 3,4 A thermocouple 
was also installed on the airfoil to verify the temperature 
on the Inconel sheet beneath the liquid crystal coating. 

Test Procedure 


low thermal conductivity material used for the airfoil 
makes interior heat conduction loss negligible, and thus 
it is neglected. The area A was the measured area of the 
Inconel heater sheet. The temperature of the test surface 
T lc was the calibrated liquid crystal temperature. The 
local recovery temperature Tj^ of the free-stream air 
around the airfoil was calculated as 

Tr = T st + r(T t - TJ (2) 

where T st and T t are the local static and total air tem- 
peratures around the airfoil, and r is the recovery factor, 
defined as the cube root of the Prandtl number. 

Reynolds Number Measurement 


The Reynolds numbers were obtained by varying 
the inlet flow rate while keeping the inlet pressure 
approximately constant. Turbulence intensity meas- 
urements were obtained by using a constant temperature 
hot film anemometer, which was installed in a traversing 
actuator on top of the test section to survey the flow 
channel in a pitchwise direction at midspan. Heat 
transfer tests were conducted for Reynolds numbers of 
0.75, 1.5, 3.0, 5.0, and 7.0xl0 6 . The inlet Mach num- 
bers ranged from 0.027 to 0.27; exit Mach numbers 
ranged from 0.058 to 0.71. 

Data for turbulence intensity, pressure distribution 
(using the pressure distribution airfoil), and heat trans- 
fer (using the heat transfer airfoil) were recorded at 
Reynolds numbers of 0.75, 1.5, 3.0, 5.0, and 7.0xl0 6 . 
The first series of tests were conducted with a clear 
tunnel (no turbulence grid installed); the second series 
used the l/8-in. grid; and the last series used the 1/4- in. 
grid. 


Data Reduction Procedure 


Heat Transfer 


Surface heat transfer coefficients were obtained 
from the energy balance 


Qi-Ql 
A(T lc - t r) 


(1) 


The Reynolds number is based on the axial chord 
and the exit conditions: 


Re = 


C(pV) exit 




( 3 ) 


where C is the chord length, p is the fluid density, V is 
the fluid velocity, and fi is the dynamic viscosity, 
(pV) exit = W/A exit , and 

= S cos 0 (4) 


S is the blade pitch, and 6 is the flow exit angle. In this 
paper, the flow angle is taken as equal to the camber 
angle (61.7°), which is approximately the case for sub- 
sonic flow. W is measured by a sharp-edged orifice 
downstream of the test section. 

Turbulence Measurements 


Turbulence intensity was measured with a TSI 
model 1210-20 hot film probe located 6.19 in. down- 
stream of the turbulence grid. Flow was normal to the 
probe axis. The probe was calibrated in the tunnel, up- 
stream of the turbulence grid and boundary layer bleeds, 
by varying the Reynolds number over the range used in 
this study. The probe voltage was plotted against the 
calculated density- velocity product, and a fourth-order, 
least-squares curve fit was obtained: 


where the local heat transfer coefficient h was calculated 
at the location of the calibrated color band (an isotherm 
which, in this case, also represents a uniform heat trans- 
fer coefficient). The heat energy Qj supplied to the 
Inconel heater sheet was calculated from the measured 
voltage across the sheet and the current through a shunt 
resistor in series with the sheet. The heat loss Q L was 
the calculated radiation loss from the surface. The very 


p v = B 0 + B,E + B 2 E 2 + B 3 E 3 + B 4 E 4 ( 5 ) 

where pV is the density- velocity product, B is the 
polynomial coefficient, and E is the mean voltage read 
from an integrating voltmeter. The turbulence intensity 
was then obtained from the calibration curve as 
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Tu = pV' = d[ P V(E)] dE (6) 

pV dE pV (E) 

where dE is taken as the root mean square of the fluctu- 
ating voltage read from a true rms meter. The turbu- 
lence intensity was measured in midchannel because it 
represented the location in the flow channel least 
affected by the airfoil leading edge and the sidewall. 

Results and Discussion 

Figure 5 presents turbulence intensities as a func- 
tion of Reynolds number for three test conditions: clear 
tunnel (no grid), 1/8- and 1/4- in. grids. In general, the 
turbulence intensities show a stronger Reynolds number 
dependency at lower Reynolds numbers. In all cases, the 
turbulence intensities decreased only slightly for 
Reynolds numbers greater than 3.0X10 6 . The turbulence 
intensities for the clear tunnel (no grid) shown here are 
higher than those obtained from a typical quiescent 
wind tunnel because of the turbulent pipe flow condition 
that exists ahead of the test section. For the cases of the 
1/8- and l/4-in. bar grids, the turbulence intensities are 
in general agreement with those calculated from Baines 
and Peterson 5 for Reynolds numbers greater than 
3.0X10 6 . 

Pressure Distribution 

Both the airfoil and the contoured sidewalls were 
instrumented with the same number of pressure taps to 
check the effect of periodicity, as well as the three- 
dimensionality, of the flow in the cascade. 

Figures 6(a) and (b) show the effect of turbulence 
intensity on the airfoil pressure distribution at midspan 
at Reynolds numbers of 3.0 and 7.0x10 , respectively. 
The turbulence levels varied from 12.3 to 1.8 percent 
depending on the Reynolds numbers and the turbulence 
grids used in these tests. The figures indicated that tur- 
bulence intensity has little or no effect on the airfoil 
pressure distribution. The turbulence intensity effects at 
other Reynolds numbers are also negligible. Similar plots 
for the airfoil hub and tip regions of the present experi- 
ment also show little or no turbulence intensity depen- 
dency on pressure distribution. These findings are in 
general agreement with Dring et al., 6 who reported on 
the aerodynamic and heat transfer characteristics of a 
large-scale rotating turbine model. They found that the 
midspan pressure distributions were essentially unaf- 
fected by the turbulence generating grid. 

Figure 7(a) compares the airfoil pressure distribu- 
tion at the hub, midspan, and tip for a Reynolds num- 
ber of 3.0xl0 6 and a turbulence intensity of 6.5 percent 
(l/8-in. grid). The figures show that, except near the 


suction surface leading edge, there is no variation in the 
spanwise pressure distribution. This is an indication that 
the flow on the airfoil is mostly two-dimensional. Simi- 
lar results are noted in Fig. 7(b) for a Reynolds number 
of 7.0xl0 6 and a turbulence intensity of 5.9 percent. 

For the l/8-in. grid case, Figs. 8(a) to (c) compare 
the pressure distribution between the airfoil and the con- 
toured sidewall at midspan for Reynolds numbers of 3.0, 
5.0, and 7.0xl0 6 . The agreement between the airfoil and 
the contoured sidewall is good on the pressure surface 
for all the Reynolds numbers shown. For the suction 
surface, the agreement is also good for about 50 percent 
of the airfoil. From about 50 percent of the chord to the 
trailing edge, there is deviation between the airfoil and 
sidewall pressure distribution. The deviation increases 
progressively with Reynolds number. The deviation is 
especially obvious at a Reynolds number of 7.0xl0 6 . It 
is speculated here that the reason for the deviation at 
the suction surface near the trailing edge is that the 
facility is not a true cascade. In a true cascade, the 
airfoil suction surface would be uncovered downstream 
of the throat. In this facility the continuation of the 
contoured sidewall pressure surface that forms the cas- 
cade wall would cause the air to continue to accelerate 
downstream of the throat, thereby causing the pressure 
to decrease. This effect is seen at all Reynolds numbers; 
however, it was especially pronounced at Re = 7.0x10 . 

Heat Transfer 

In Fig. 9(a) the heat transfer coefficient is plotted 
against the normalized surface distance s/L for Reynolds 
numbers of 0.75, 1.5, 3.0, 5.0, and 7.0xl0 6 for the clear 
tunnel (no grid) case. As seen in the figure, heat transfer 
increases with Reynolds number, as expected. On the 
pressure surface, the transition (as evidenced by the 
sharp rise in the heat transfer coefficient) appears to 
occur near an s/L of 0.09, which is near the point of 
tangency of the leading edge circle. Examination of the 
liquid crystal photographs reveals the presence of a weak 
separation at this point. The pressure distribution con- 
tour also indicates a small adverse gradient near this 
region. It is speculated that flow separation at this point 
triggers the transition for the three highest Reynolds 
numbers; this is the reason why the transition point does 
not change with increasing Reynolds number. After the 
initial rise due to transition, the heat transfer coefficient 
continues to increase at a reduced rate; this is caused by 
the flow acceleration on the pressure surface. 

On the suction surface, the location of the transi- 
tion moves upstream with increasing Reynolds numbers, 
indicating that no flow separation is present. This is 
supported by the fact that the pressure distribution on 
the suction surface does not show any adverse pressure 
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gradient and that the liquid crystal photographs do not 
show evidence of flow separation. 

Near the suction surface trailing edge, the heat 
transfer coefficient for a Reynolds number of 7.0 xlO 6 
shows an increasing trend, which is incompatible with 
the diffusion in this region shown in Fig. 8(c). The 
reason for this anomaly is not known at this time. 

Figures 9 (b) and (c) present the heat transfer 
coefficients for the moderate and high turbulence cases 
(1/8- and 1/4- in. grids). On the pressure surface, the 
separation at s/L = 0.09 is still present and can be seen 
as a small dip in the heat transfer coefficient at the two 
lowest Reynolds numbers. As the Reynolds number is 
increased, the transition moves upstream from this 
point. For the two highest Reynolds numbers, there is 
an abnormal drop in the heat transfer coefficient at the 
end of the transition region that corresponds to the 
separation point. This drop was also observed by 
Hippensteele, Russell, and Torres. 4 

In comparison with Fig. 9(a) for the clear tunnel 
case, heat transfer for the suction surface and at the 
leading edge region is higher because of the higher tur- 
bulence intensity resulting from the 1/8- and l/4-in. 
grids used. Also, transition from laminar to turbulent 
flow occurs sooner and over shorter distances. For the 
same reason, heat transfer for both the pressure and 
suction surfaces is also higher than that for the clear 
tunnel case. These trends were also noted by Dring 
et al. 6 

At the leading edge, where the flow is laminar, 
heat transfer varies as Re 0,5 . A Frossling number (Nu d / 
(R e d inlet) 0 * 5 ) plot is useful: Frossling numbers for 
similar turbulence levels should fall on the same curve. 
Figure 10(a) is a plot of the Frossling numbers at the 
leading edge region for the clear tunnel case. The figure 
shows that, with some exceptions, the Frossling numbers 
in the leading edge region do generally fall on a single 
curve. The reason why some Frossling numbers are 
lower is not known at this time. Figure 10(a) also shows 
that there is a sharp increase in heat transfer in this 
transition regime, especially for the higher Reynolds 
numbers. Aft of the leading edge region, heat transfer 
departs from laminar behavior and will no longer vary 
as the square root of the Reynolds number. 

Figure 10(b) shows Frossling numbers with moder- 
ate turbulence intensity (l/8-in. grid) upstream of the 
vane. Frossling numbers for the leading edge region 
again generally fall on a single curve. Because the stag- 
nation region heat transfer is augmented by higher free- 
stream turbulence intensity, the level of this curve is 
higher than that for the lower turbulence case of 
Fig. 10(a). Figure 10(b) also shows that transition from 


laminar to turbulent flow also occurs sooner, especially 
for those at higher Reynolds numbers. Figure 10(c) 
shows Frossling numbers for high turbulence intensity 
(l/4-in. grid) upstream of the vane. The characteristics 
are similar to those of Fig. 10(b). 

It would be reasonable to raise the question that, 
if laminar heat transfer at the leading edge with similar 
turbulence intensity levels can be normalized by the 
Frossling number, is it possible, for similar turbulence 
intensities, to normalize the turbulent heat transfer for 
the remaining airfoil surfaces? One possible approach 
would be to use the Nusselt number (based on the dis- 
tance from the stagnation point) divided by the 
Reynolds number (based on the axial chord and exit 
conditions) raised to the 0.8 power (Nu g /(Re c exit ) 0 * 8 ). 
Such an attempt was made using data from the 1/8-in. 
grid (see Fig. 11). In this figure, heat transfer values at 
the leading edge should be disregarded because the 
plotting parameter is not valid for the leading edge area. 
For the remaining airfoil surfaces, the figure shows a 
data spread of approximately ±6 percent on both the 
suction and pressure surface trailing edges. The data 
spread is larger near the transition region because 
transition occurs at different locations and over different 
transition lengths at different Reynolds numbers. 

From the foregoing discussions, it may be con- 
cluded that, for similar turbulence intensities, heat 
transfer for high Reynolds numbers at the leading edge 
can be correlated by using the Frossling equation devel- 
oped for lower Reynolds numbers. For the rest of the 
airfoil, heat transfer for similar turbulence intensities 
may be approximately correlated by using the parameter 
Nu s /(Re c , exit ) 0 - 8 . 

Figure 12(a) compares the effect of turbulence in- 
tensity at Re numbers of 0.75 XlO 6 for the clear tunnel, 
1/8- and 1/4- in. grid cases, respectively. The increase in 
stagnation region heat transfer from the clear tunnel to 
the 1/8-in. grid (higher levels of turbulence) case is very 
evident. However, there seemed to be no difference be- 
tween the 1/8- and 1/4- in. grid cases at the leading edge 
region. The fact that the heat transfer coefficient for the 
1/8- and the 1/4- in. grid cases falls on the same curve is 
somewhat unexpected, considering the l/8-in. grid pro- 
duces a turbulence level of about 6 percent and the 
1/4- in. grid about 11 percent. The reason for this result 
is not known, but it could possibly be attributed to the 
larger turbulent length scale for the 1/4- in. grid. The 
heat transfer coefficient in the stagnation region is pro- 
portional to the turbulence intensity, and it has been 
shown to be inversely proportional to the ratio of the 
length scale to the leading edge diameter (G. James 
VanFossen and Chan Y. Ching, “Measurement of the 
Influence of Turbulence Integral Length Scales on Stag- 
nation Region Heat Transfer,” to be published in 1993, 
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NASA Lewis Research Center, Cleveland, Ohio). The 
figure also shows that higher turbulence levels move 
transition further upstream. 

Figure 12(b) compares the effect of turbulence 
intensity at Re = 7.0xl0 6 . At this high Reynolds num- 
ber, the transition point has moved upstream almost to 
the stagnation point. The heat transfer coefficients for 
the 1/8-in. grid are still higher than those for the clear 
tunnel (no grid), but the heat transfer coefficients for 
the 1/8- and the 1/4- in. grid cases are no longer coinci- 
dent in the stagnation region; the l/4-in. grid has the 
highest heat transfer coefficients there. 

In Fig. 13 the Frossling number at the stagnation 
point is plotted against a parameter developed by 
VanFossen and Ching. This parameter involves turbu- 
lence intensity, Reynolds number, and integral length 
scale. The solid curve is their correlation and the dashed 
curve is a ±8-percent error band. Length scales for the 
present data were estimated from VanFossen and 
Ching’s curve fit of the length scale versus the distance; 
they used 1/2-, 1/4-, 1/8-, and 1/16-in. grids which had 
the same ratio of bar width to mesh spacing as that of 
the present test. For the data with clear tunnels (no 
grid), the length scale was set to the width of the test 
section, 8.5 in. Agreement with the correlation is reason- 
able especially at the high Reynolds numbers for the two 
turbulence grids. Agreement for the clear tunnel case is 
not quite as good but is still reasonable considering the 
greater uncertainty in the length scale used. 

Figure 14 compares the pressure distributions of 
the airfoil and the contoured sidewall at Re = 7.0 xlO 6 
with the computed pressure distribution from 
TRAFC2D superimposed. The predicted values agree 
well with the contoured sidewall experimental data, ex- 
cept that near the suction surface leading edge area. The 
small disagreement on the suction surface leading edge 
area is probably caused by the three-dimensional effect. 
Because the code used was a two-dimensional version, it 
was not able to follow the three-dimensional pressure 
distribution. The agreement is good for the rest of the 
sidewall surface. Agreement is generally good with the 
experimental pressure distribution on the airfoil, except 
in certain areas: near the suction surface leading edge, 
there is again the possible three-dimensional flow field 
effect, as just discussed; near the suction surface trailing 
edge area, the airfoil pressure was decreased because of 
the presence of the contoured sidewall (see the earlier 
discussion in the Pressure Distribution section). 

A comparison of the predicted and experimental 
heat transfer coefficients is presented in Fig. 15 for 
Reynolds numbers of 1.5, 3.0, and 7.0xl0 6 for the clear 
tunnel case. The predictive tool used was again 


TRAFC2D. As can be seen from the figure, the heat 
transfer coefficients on the suction surface are relatively 
well predicted. On the pressure surface, the heat transfer 
coefficients are uiiderpredicted for all the Reynolds num- 
bers. The predicted heat transfer values show encourag- 
ing results; however, it is evident that additional work 
on the predictive techniques is required. 

Uncertainty Analysis 

An uncertainty analysis, based on the work of 
Abernathy and Benedict 7 showed that the uncertainty 
is a function of Reynolds numbers. For all Reynolds 
numbers and all turbulence levels tested in this facility, 
a conservative estimate of the uncertainty of the heat 
transfer coefficients is approximately 6 percent, with the 
exception of two cases (Re = 3.0 and 5.0xl0 6 , for the 
clear tunnel). The uncertainty of the heat transfer coef- 
ficients for these two worst cases is 11 and 8.7 percent, 
respectively. The high uncertainty level of these two 
cases is attributed to the extremely warm air tempera- 
ture in these tests (91 °F), resulting in a temperature 
difference of only 10 °F between the air and the liquid 
crystal. 

Conclusions 

Aerodynamic and heat transfer data were presented 
for Reynolds numbers ranging from 0.75 to 7.0xl0 6 and 
turbulence intensities ranging from 1.8 to 15.1 percent. 
The range of Reynolds numbers encompassed those en- 
countered in current aeropropulsion turbines and in the 
Space Shuttle main engine (SSME) turbopump turbines. 

From the pressure distribution plots presented, it 
can be concluded that 

1. The flow in the cascade was nonperiodic at the 
airfoil suction surface near the trailing edge. This 
nonperiodicity was evident in all Reynolds numbers 
although it was not very pronounced except at Re 
= 7.0xl0 6 . 

2. The flow over the airfoil was primarily two- 
dimensional. 

3. The turbulence intensity had a minimal effect on 
the pressure distribution, as expected. 

4. A numerical solution using the TRAFC2D code 
for pressure distribution agreed well with the contoured 
sidewall data. 

With respect to heat transfer, the results of this 
study indicated that 
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1. Reynolds number had a large effect on heat 
transfer and on the transition from laminar to turbulent 
flow. 

2. For the high Reynolds numbers typical of the 
Space Shuttle main engine turbopump, heat transfer was 
dominated by transitional and turbulent flow down- 
stream of the stagnation region. 

3. For a given turbulence intensity level, heat 
transfer at the leading edge for all Reynolds numbers 
can be correlated with the Frossling number developed 
for lower Reynolds numbers. 

4. For a given turbulence intensity, heat transfer 
on the airfoil downstream of the leading edge can be 
approximately correlated by using a dimensionless 
parameter (Nusselt number, based on the distance from 
the stagnation point, divided by the Reynolds number, 
based on the axial chord and exit conditions, raised to 
the 0.8 power) to within ±6 percent at the trailing edge. 

5. A comparison of the experimental heat transfer 
data with the numerical solution of the TRAFC2D code 
showed encouraging results. 
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Figure 1 . — Variable Reynolds Number Heat Transfer Facility. 



Figure 2. — Test section. 



C-88-11421 

Figure 3.— Pressure distribution airfoil showing pressure tap 
locations on suction surface. 
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Figure 4. — Heat transfer airfoil showing grid markings for the 
locations of temperature distribution. 



Figure 5. — Turbulence intensity as function of Reynolds number 
for three test conditions. 
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Figure 6. — Effect of turbulence Intensity on airfoil pressure 
distribution at midspan. 
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Figure 7. — Airfoil pressure distribution at midspan, hub and 
tip. 
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Figure 8. — Comparison of airfoil and contoured sidewall 
pressure distribution at midspan fori /8-in. grid. 
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Figure 9. — Effect of Reynolds number on heat transfer coefficient. 
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Figure 1 0. — Frossling correlation of variable Reynolds number heat transfer data at leading edge. 
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Figure 1 1 .—Correlation of variable Reynolds number heat transfer data for 1/8-in. grid case. 
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Figure 14. — Comparison of computed and experimental 
pressure distributions of airfoil and contoured sidewall at 
Reynolds number of 7.0x1 0 6 . 
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Figure 1 5. — Comparison of computed and experimental airfoil heat transfer coefficients for clear 
tunnel (no grid) condition. 
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